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Cold neutral atoms provide a versatile and controllable platform for emulating various quantum systems. 
Despite efforts to develop artificial gauge fields in these systems, realizing a unique ideal-solenoid-shaped 
magnetic field within the quantum domain in any real-world physical system remains elusive. Here we 
propose a scheme to generate a "hairline" solenoid with an extremely small size around 1 micrometer which 
is smaller than the typical coherence length in cold atoms. Correspondingly, interference effects will play a 
role in transport. Despite the small size, the magnetic flux imposed on the atoms is very large thanks to the 
very strong field generated inside the solenoid. By arranging different sets of Laguerre- Gauss (LG) lasers, the 
generation of Abelian and non-Abelian SU(2) lattice gauge fields is proposed for neutral atoms in ring- and 
square-shaped optical lattices. As an application, interference patterns of the magnetic type-I 
Aharonov-Bohm (AB) effect are obtained by evolving atoms along a circle over several tens of lattice cells. 
During the evolution, the quantum coherence is maintained and the atoms are exposed to a large magnetic 
flux. The scheme requires only standard optical access, and is robust to weak particle interactions. 

There is intense interest in finding new settings in which different forms of gauge fields appear. Ultracold 
atoms in optical lattices are ideal systems with which to achieve this goal as they offer unprecedented 
possibilities of emulating condensed matter and high energy physics models^"^^. The original magnetic 
and electric AB effects (type-I) are distinguished from the neutral scalar AB effect and the Aharonov-Casher 
effect (type-II) by the absence of any electromagnetic fields^^. This original AB effect shows that electron wave 
packets are influenced by the non-zero potentials and can obtain a non-zero phase shift although they travel in 
field-free regions^^. The type-I AB effect is important conceptually because it bears on three issues: whether 
potentials are physical quantities, whether action principles are fundamental, and the principle of locality. As such 
it has spawned many studies^^"^^. In 1986, to overcome the issue of stray fields in the previous experiments, 
Toromura et al. have performed a beautiful experiment to prove the presence of the AB phase by using magnetic 
toroids with superconducting shields to eliminate the leakage fields^^ In 2007, Adam et al. tested the magnetic 
type-I AB effect in the absence of any force for a macroscopic system^^. Although the phase shift has been 
observed experimentally, the direct interference signatures in the original version of the AB effect have not yet 
been reported. In solid-state AB interferometer, the AB effect has also been nicely demonstrated^^ Given the great 
importance of type-I AB effect as a quantum mechanical phenomenon, it is of fundamental interest to ask 
whether the original AB effect can be reproduced in a quantum regime and how to observe the associated 
interference effects. A good realization in degenerate cold gases has however so far proved elusive. 

For the simulation of uniform or staggered gauge fields, several studies have proposed a number of methods^"^, 
including rotation of the confining trap^°"^^, using adiabatic motion of multilevel atoms in a dark-state picture^^"^^ 
and employing spatially dependent two-photon dressing^^'^^. In the presence of a lattice potential, schemes to 
generate gauge fields have been proposed using staggered- aligned states^^"^^ or through rotations^^. In this work, 
we propose a scheme to emulate ideal-solenoid-shaped Abelian and non-Abelian gauge fields by employing LG 
lasers^^. The generated magnetic field is strongly localized within a lattice cell, which allows an initially prepared 
coherent state of atoms to evolve along a circle without any electromagnetic field. Due to the small radius of the 
solenoid, the atoms can evolve along two paths within coherence time and eventually overlap so that interference 
occurs. At the same time, the strong strength of the magnetic field exposes the atoms to a large magnetic flux 
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which results in a clearly visible interference signature. Besides mir- 
roring the original AB effect relating to an Abelian gauge field, we 
also investigate systems exposed to a non-Abelian field which is 
achievable by employing more lasers. Different from the AB phase 
mentioned in other proposals which is accumulated from atomic 
evolution along four sides of a lattice cell, the AB phase here is 
accumulated from atomic evolution along a circle enclosing several 
tens of lattice cells, which makes the detection of interference more 
readily accessible to experiments. 

Results 

In this work, we consider a scheme to realize a "hairline" solenoid in 
both a square lattice [Fig. 1(a)] and in a ring geometry^^ [Fig. 1(b)] by 
employing LG-laser-assisted tunnelling. Our scheme consists of 
three ingredients. Firstly, atoms in different internal states are 
trapped in a staggered manner in a ring or a square lattice, e.g., as 
shown in Fig. 1(c) for a square -lattice case. The open and solid circles 
represent atoms in different internal states, respectively. Secondly, 
the depth of lattice potentials is tuned to be very large such that 
regular tunnelling among lattice sites is prohibited. Finally, the addi- 
tional lasers with non-zero orbital angular momentum are switched 
on to induce tunnelling between the adjacent lattice sites. 

In particular, we consider a trapping scheme with the basic prin- 
ciple proposed in previous works^^'^°, where two-electron atoms gen- 
erally possess a spin -singlet ground state g and a long-lived spin 
triplet excited state e with a lifetime around 20 s^^'^l The g and e 
states have opposite polarizabilities at the "anti-magic" wavelength 
X^l ^m^°, where they experience potentials with opposite signs and 
therefore align in a staggered manner, as shown in Fig. 1(c) for the 
square-lattice case. We would like to note that apart from the typical 
square lattice system which can be created by the interference of 
counter-propagating laser beams at the "anti-magic" wavelength, a 
ring lattice in (b) can also be created by interfering an off- resonant 
LG laser and a plane wave with their wavelengths chosen to be at the 
"anti-magic" wavelength^l At this stage, the lattice depth is set to be 
very large such that the tunnelling of atoms is strongly suppressed. As 
shown in Fig. 1(d), resonant LG lasers are appUed to drive transitions 
between the g and e states, leading to atomic hoppings along a loop as 
shown in Fig. 1(c). Different values of / correspond to different accu- 
mulated phase when atoms move around the loop. We consider / = 1 




Figure 1 | Schematic diagram for generating strongly localized effective 
Abelian and non-Abelian gauge fields with cold atoms trapped in square 
(a) and ring (b) lattices. The atoms can hop from one site to a neighbouring 
site via laser-assisted tunnelling by shinning LG lasers perpendicular to the 
lattice surface (a). The LG beams drive the spin-conserved transitions 
between two identical sublevels in neighbouring sites in (d) to generate 
Abelian fields, and spin-flipping transitions between two different 
sublevels in neighbouring sites in (e) to generate SU(2) non-Abelian fields. 
The resulting magnetic flux as shown in (c) is non-zero when going around 
a loop (e.g., the blue square) including the laser centre which is addressed 
in the lattice centre and zero when the laser centre is excluded from the 
loop (e.g., the green square). 



and / = 0 lasers in this work. For the former with / = 1, it is found that 
the accumulated phase is n, which resembles a system where an 
Abelian magnetic field is applied to the atoms. For the latter with 
/ = 0, the accumulated phase is 0, thus emulating a system with a zero 
magnetic field applied to the atoms. Moreover, when considering two 
sub-states in each ^ and e levels, which is shown in Fig. 1(e), an SU(2) 
field can be generated by employing two LG lasers with / = 0 and 
/ = 1 to drive transitions between different sub-states. 

We first consider a Hamiltonian describing atoms with only one 
sub -state in each g and e levels as 



H=J2 [dr^li 

s=g,e J 



11 

2m 



dt\d,,E{rye""'^i,\{r:)^|tJr) + h.c. 



gs,s' y'rfr^J(r)^:,(r)^,(r)^,(r) 



(1) 



The four terms in turn give the kinetic energy, the lattice potential, 
the laser-assisted transition, and the interaction between atoms, 
respectively. Here, ij/sir) represents the atomic field operator at posi- 
tion r, p is the momentum operator, and m is the mass. The state- 
dependent sign of the lattice potential y (r) is denoted by rjg= + and 
r]e= ~, due to the lattice lasers at the "anti-magic" wavelength. The 
dipole moment of the g-e transition is labefled as deg, and the inter- 
action strength is given by ^5^5'. As shown in Fig. 1(a), the LG laser is 
applied perpendicular to the lattice. The amplitude of the LG laser 
resonant to the g-e transition reads 



E{r)=Efpi{r)e^^^e^ 



'i(p i{(Dt — kz) 



where /p,(r) = (-l) 



2p\ 



7i(p+|/|)! 



x/2r 



(2) 



, is the 



waist of the beam, and ijj' are the Laguerre functions^^. For an Ns 
X Ns 2D lattice, we choose = Nsa/2. The cylindrical coordinate (r, 
(p, z) is chosen that the longitudinal axis z is along the propagation 
direction of the LG laser, and the labels p and / represent the radial 
and azimuthal quantum numbers, respectively. 

We choose a lattice potential V that is minimized at sites 
0 = {open-circle sites } and maximized at sites £^ = {solid-circle 
sites} (see Fig. 1). In the presence of a deep lattice potential, xjjsi^) 
can be expressed by Wannier functions in the lowest band as 

'A^(e)W- E;G^(o4^*(^~^j)' ^^^^^ ~ Wannier 
function^^. In the numerical simulation, we approximate the 
Wannier functions to be Gaussion functions. Here, we have assumed 
that the lattice potential is symmetric for the g and e states. The 
Hamfltonian then reduces to a tight-binding model as 



-Jo ^ (^a]aj + h.c. 



(3) 



+ ^eia]ai+j^a]a]ar 



Here, the chemical potential is not relevant in terms of the gauge 
fields, so we treat it as uniform by tuning the lattice potential. We 
assume the on-site interaction strength to be identical for each site 
and label it as L^. For the first hopping term, (/, j) denotes two nearest- 
neighbouring sites that one is in the set Q whfle the other is in the set 
8, between which the tunnelling is induced by the LG laser and has a 
strength^^'^^ 
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Jij dggE 



' J drw*{r- 



■^i)fpi{r 



(4) 



where / e Q and; e 8. For the second hopping term, ((/,;)) denotes 
two next-nearest-neighbouring sites, which are in the same set Q or Sy 

and the corresponding tunnelUng strength Jq = — j drof {r — Xi) 

[p^/2m + f/5y(r)](L)(r — Tj). By controlUng the lattice potential, the 
next-nearest-neighbouring tunnelling Jq can be tuned to be much 
smaller than the nearest-neighbouring // which allows us to neglect 
/o and only focus on Ji^j in the following. 

A natural property of the laser with a non-zero orbital angular 
momentum is a spatial-dependent phase. Since the laser is appUed 
perpendicular to the lattice, the phase of Ji^j will only depend on the 
azimuthal angle cp, as shown in Eq. (4). The phase of //j can be 
estimated as //^ ^ e"''^ with cp the azimuthal angle of the midpoint 
between two adjacent sites / e Q and j e 8. We would like to remark 
that the sign of the phase depends on the tunnelling direction, where 
for / G 8 and j e Q, Jij cxd e'^^. Due to the vortex of the LG laser, the 
accumulated phase for atoms moving around a loop enclosing or 
excluding the laser centre are different. As a simplified illustration, 
we consider a centre cell of the square lattice [blue square in the 
centre of Fig. 1 (c)]. Here, we assume that the laser centre coincides 
with the centre of this cell. As a four-site circle, G and 8 sites appear 
alternatively. Without the loss of generality, we assume cp = 0 along 
the upward direction in Fig. 1 (c). Then, for a particle moving in a 
counterclockwise direction, it undergoes tunnelling with phases 0, 

y, — /tt, and ^ in the four links, respectively. As a result, the 

accumulated phase around the centre cell is given by 

OH ln-\ = In. By choosing an odd /, the accumulated phase 

is nontrivial, giving a non-zero gauge field within the centre cell. In 
the following, we will focus on the case with / = 1 for the generation 
of an Abelian gauge field. As a comparison, a laser with / = 0 will give 
a zero accumulated phase and therefore correspond to a zero gauge 
field. On the other hand, for a cell far away from the laser centre 
[green square in Fig. 1 (c)], the accumulated phase is always zero for 
any /, which means that a non-zero gauge field is strongly localized 
inside the centre cell for / = 1, resembling a very thin solenoid. As an 
illustration, we consider a cell centred around Tq. The laser centre is 
set to be the origin of the system. The coordinates of each site belong- 
ing to the cell is r = Tq + ^r, which gives the arimuthal angle 



) = arctan 



cos^o^y ~ sin^o^-^ 



ro 



provided Vq > a. Here cpo is the azimuthal angle of Tq, and a is the 

a 

lattice constant. With Sx, 3y= ± - , the accumulated phase is zero 

for a cell away from the laser centre. 

Numerical simulation for the accumulated flux (/>y is shown in 
Fig. 2, where an / = 1 LG laser drives the tunnelling of atoms along 
a loop enclosing the lattice centre that coincides with the laser centre. 
The flux is defined as = arg(^^) - arg(/,y+i) + arg(/,+i,^+i) - 
arg(//+ij), where the tunneling strength //^ is given in Eq. (4). As 
shown in Fig. 2, the gauge field is non-zero within only a few cells 
around the centre. We would like to remark that even if the centre of 
the laser is slightly shifted from that of the lattice, it will not affect the 
results significantly, although some small fluctuations may appear 
around the solenoid. In a similar way, an LG laser drives tunnelling of 
atoms which are trapped in a ring lattice. For an / = 1 laser, the 
accumulated phase along the ring is n. It is straightforward to extend 
to the SU(2) case, where ij/s becomes a 2 X 1 column matrix to 
include the two sub-states and g2 {ei and 62) in g {e) level [see 
Fig. 1 (e)] . Two laser beams with / = 0 and / = 1 are employed to drive 
g\-e2 and ^2-^1 transitions, respectively. To give a unitary hopping 




Figure 2 | Numerically calculated net flux for each cell on a 2D square 
lattice. We have chosen / = 1 for the LG laser, and the laser centre is 
addressed in the lattice centre. The flux is n when atoms tunnel along a loop 
including the laser centre. Although small fluctuations appear around the 
centre cell, the positive and negative phases cancel each other and the total 
flux for a loop including the centre cell remains a value of %. 



matrix for two spins, we should choose / = 0 and / = 1 lasers with the 
same amplitude to drive the spin-flipping transitions. It is satisfied 
when we apply an LG mode with/> = 0, / = 1, and a superposition of 
two LG modes withp = 0, / = 0 andp = 1, / = 0. The non- Abelian 
tunneling matrix is then given by 



7.7 = 



0 \k 



hi 



0 



(6) 



where Jy is given in Eq. (4). 

We would like to note the difference between our scheme and the 
works for continuous systems investigated in Refs. [13,14]. Although 
both employ laser beams with orbital angular momentum, the latter 
schemes deal with dark states in the electromagnetically induced 
transparency configuration, where the adiabatic motion of multilevel 
atoms is appUed. In the present scheme, the necessary condition for 
the generation of gauge fields is a staggered distribution of different 
internal states and laser- assisted tunnefling. Therefore, unlike a uni- 
form field generated in Refs. [13,14], a solenoidal field is produced 
here. 

Type-I AB effect. To see the interference patterns, we load a Bose- 
Einstein condensate (BEC) initially away from the centre of the LG 
laser, as seen in the plot shown in Fig. 3(a)-(c) at t = 0 for a ring- 
lattice case and in Fig. 4(a) for a square-lattice case. The system is 
then aflowed to evolve under zero, Abelian and non-Abelian gauge 
fields. And finally the fringe patterns are measured with time. After 
the initial BEC state is prepared, the LG beams are switched on. At 
this stage, the Hamiltonian governing the time evolution is 



H = - ^ Ijijolaj + h.c.]. 

(Uj) 



i7) 



Here we have neglected the interactions between atoms. Evolutions 
for weakly interacting gas are analyzed in the Methods section, where 
the results show that the key features of the interference patterns 
remain robust against a weak interaction. For an LG laser with / = 
1, the accumulated phase along the circle is n. Since the atoms move 
along left or right path, the atoms from different paths wifl possess a 
different phase factor at the opposite site, giving a destructive 
interference. For an LG laser with / = 0, the phase is always zero, 
which resembles the system with particles moving in the absence of 
any gauge field. The atoms evolving along two different paths wifl 
possess the same phase factor at the opposite site, and the 
interference is destructive. Therefore, the interference fringes are 
closely related to the phase in our scheme. 
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Figure 3 | The time evolution of particle distributions for particles 
hopping around a loop formed by LG beams in a ring with ATl = 100 sites 
under a zero gauge field (a), an Abelian U(l) field (b), or a non-Abelian 
SU(2) field ((c) and (d)). The particle densities are represented by the 
percentages of particle numbers to the total number of the initially 
prepared Bose-Einstein condensate. At time t = 0, the particles are 
prepared around the site x = 1. During the time period from t = 0 to 
t = T = 600h/ERy the fringe patterns around the opposite site x = NJ2 
become completely different for the three field cases (B = 0, U(l) and 
SU(2)). Around this site, there is always a constructive interference for the 
/ = 0 case in (a) and a destructive interference for / = 1 in (b). In contrast, 
for the non-Abelian case, the charge wave density, the sum of the two 
effective spins, has no unique behavior as shown in (c), and the spin wave 
density (d), as a difference of the two effective spins, exhibits the red-color 
(blue-color) component with a constructive (destructive) interference 
around the site x = NJ2 and a destructive (constructive) interference 
around the nearest neighboring area. Movies of the time evolution are 
available in the supplementary materials. 

For a ID ring-shaped optical lattice as illustrated in Fig. 1(b), 
numerical simulations for spatial distributions of particle numbers 
are plotted in Fig. 3, where the particle numbers are shown as a 
percentage of the total number in the initially prepared BEC. We 
have chosen the lattice site number as Ni^= 100 and the time period 
T as 600 in unit of Er/H, where Er is the recoil energy and the typical 
hopping is Jc^O.OSEr. With ER/n = 2n x 900 Hz^\ the time scale is 
approximately in the range of 100 ms. Distinctly different interfer- 
ence pattens can be seen around x = NJ2 site for particles experi- 
encing a zero magnetic field as shown in (a), an Abelian U(l) gauge 
field as shown in (b), and a non-Abelian SU(2) gauge field as shown 
in (c)-(d). The reason for the two subfigures in the SU(2) case is that 
two effective spins are involved in the SU(2) field. We plot an effec- 
tive charge wave density (the sum of densities of the two effective 
spins) in (c) and an effective spin wave density (the difference of 
densities of the two effective spins) in (d), respectively. At the NJ2 
site, a destructive (constructive) interference is always seen in (b) 
[(a)]. Two components appear with one colored red and the other 
blue, where the constructive interference occurs for one type (red). 

For the 2D square lattice case as illustrated in Fig. 1(a) and 1(c), we 
can also prepare the initial state by loading a BEC around one pla- 
quette. For the U(l) Abehan field case as shown in Fig. 4(b), we turn 
on an / = 1 laser as illustrated in Fig. 1(d), and for the non-Abelian 
SU(2) field case as shown in Fig. 4 (c)-(d), we turn on two ortho- 
gonally polarized / = 0 and / = 1 lasers as illustrated in Fig. 1(e). 
Numerical simulations with time show that the angular momentum 
imparted by the LG beams results in a circular distributions of atoms. 
The particle number distribution at the opposite side of the circle 




Figure 4 | The time evolution of particle distributions for particles 
hopping around a loop formed by LG beams in a 2D square lattice with 40 
X 40 sites. At time t = 0 as shown in (a), the particles are distributed 
around one site. The screenshots at time t = 4000?z / Er are given for the 
system under an Abelian U(l) (b) or a non-Abelian SU(2) ((c) and (d)) 
gauge field. The particle densities are represented by the percentages of 
particle numbers to the total number of the initially prepared Bose- 
Einstein condensate. While a destructive interference always occurs 
around the opposite site in the Abelian field case (b), a nonvanishing 
charge wave density occurs everywhere in (c) and spin waves appear in (d) 
with red- and blue-colored components, where the red-color (blue-color) 
component has a constructive (destructive) interference at the opposite 
side and a destructive (constructive) interference at its two neighboring 
area. Movies of the time evolution are available in the supplementary 
materials. 



from the initially prepared site again shows clear destructive inter- 
ference in the case of the Abelian gauge field - the signature of the AB 
effect. As a concrete example, screenshots at time t = 4000h/ER are 
given in Fig. 4 (b)-(d). The effect is similar to the ring geometry case. 

Our key result for the ring geometry is a clear demonstration of the 
magnetic type- 1 AB effect, where in the presence of an Abelian gauge 
field, the AB phase always induces destructive interference at the NJ 
2 site. Absence of any population at the NJ2 site at any time is 
therefore a clear signature of the Abelian AB effect. Furthermore, 
we demonstrate the effects of the non-Abelian gauge field in terms of 
the spin density waves. 

Discussion 

We have presented a method to generate a "hairline" solenoid with 
an extremely small size around 1 micrometer. Clear evidences of AB 
effect in the interference patterns for both Abelian and non-Abelian 
SU(2) fields are obtained with a unique topological structure in both 
a ring geometry and a square lattice. A maximal magnetic flux is 
achieved as well as quantum coherence maintained during the evolu- 
tion process of atoms due to the strong strength of generated mag- 
netic field and the small radius of the solenoid. In the detection area, 
the interference is always constructive for the case of zero field and 
always destructive for the Abelian field case. For the non-Abelian 
field case, the spin wave density exhibits two components with a 
constructive or destructive interference in the detection area. The 
scheme requires standard optical access within the reach of current 
techniques, and the key conclusions are not affected by the weak 
particle interactions. Moreover, these effects are robust against deco- 
herence as they are topological in origin. Consequently, they could in 
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principle be harnessed for quantum information processings^ The 
technique is also straightforwardly generalized to SU{N) gauge fields. 

Methods 

Interaction effect. In the main text, we have considered the non-interacting gas. A 
weakly interacting degenerate Bose gas can be described in the framework of the 
Gross -Pitaevskii equation, and now we will treat our models at the Gross-Pitaevskii 
level. To start with, we write down the Hamiltonian that describes interacting bosons 
evolving on a lattice and subjected to a gauge field: 

H = - ^ (ji^ja]aj-\-h.c.] + — ^a]a]aiai. (8) 

Here a, is a canonical Bose annihilation operator on sites of the optical lattice labeled 
by the integer z, is the hopping amplitude, and U is the interaction strength. For a 
weakly interacting gas, which is the typical case for the atomic gas'^'', the low 
temperature dynamics of H can be described by introducing the complex dynamical 

mean field variable whose value is a measure of j \/Jri). Here (n) is the 

average atom occupancy per lattice site. The dynamics is then described by the 
discrete Gross-Pitaevskii equation^" 



where /I = . By numerically solving the differential equations for evolutions 

driven by systems exposed to zero, Abelian and non-Abelian gauge fields, we have 
shown that the key features of the interference patterns presented here remain robust 
up to at least /I = 5, which is easily achievable within current experimental setups. 
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